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ABSTRACT 



The bulk of recent cosmological research has focused on the adiabatic cold dark matter model 
and its simple extensions. Here we present an accurate fitting formula that describes the matter 
transfer functions of all common variants, including mixed dark matter models. The result is a 
function of wavenumber, time, and six cosmological parameters: the massive neutrino density, 
number of neutrino species degenerate in mass, baryon density, Hubble constant, cosmological 
constant, and spatial curvature. We show how observational constraints — e.g. the shape of the 
power spectrum, the abundance of clusters and damped Lya systems, and the properties of the 
Lya forest — can be extended to a wide range of cosmologies, including variations in the neutrino 
and baryon fractions in both high-density and low-density universes. 



Subject headings: cosmology: theory - dark matter - large-scale structure of the universe 



Introduction 



Most popular cosmologies rely on density perturbations generated in the early universe and amplified 
by gravity to produce the structure observed in the universe, such as galaxies, galaxy clustering, and the 
anisotropy of the microwave background. It is of particular interest that the spectrum and evolution of 
these fluctuations depends on the nature of the dark matter as well as upon the "classical" cosmological 
parameters. Hence, by the study of the observable signatures of the perturbations, one hopes to learn not 
only about quantities such as the density of the universe or the Hubble constant, but also what fraction of 
the matter in the universe is in baryons, cold dark matter (CDM), massive neutrinos, and so forth. 

The calculation of how the various types of dark matter and the background cosmology affect the power 
spectrum can be treated for much of the history of the universe using linear perturbation theory. Numerically, 
this reduces to integrating the coupled Boltzmann equations for each mode as a function of time. For modes 
above the Jeans scales of the gravitating species, the growth of perturbations is independent of scale. In 
the absence of hot or warm dark matter, this scale drops precipitously after recombination and therefore 
the late-time power spectrum may be separated into a function of scale and a scale-independent growth 
function that incorporates the effects of time, cosmological constant, and curvature. These growth functions 
are well-cataloged (e.g. Peebles] 1980), while the form of the spatial function can be found numerically (e.g. 
Bond fc Efstathiou 1984) or quoted from fitting forms (e.g. Bond fc Efstathiou 1984; Bardeen et al. 1986; 
Holtzmarj 1989; Eisenstein & Hu 1997, hereafter [EH971) . 

With the introduction of massive neutrinos (Fang et al. 1984; Valdarnini fc Bonomctto] 1985; Achilli et 
al. 1985), or other forms of hot dark matter, the spatial and temporal behavior of the perturbations cannot 
be separated. The Jeans scale, also called the free-streaming scale, of the neutrinos remains significant 
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after recombination (Bond & Szalay 1983). In this case, the CDM and baryon perturbations are not traced 



by neutrinos and grow more slowly due to the reduction of the gravitational source (Bond et al. 1980). 
Since the free-streaming scale itself moves to ever smaller scales with increasing time, the transfer function 
acquires a non-trivial time dependence. Consequently, a cosmological constant or non-zero curvature enters 



the problem in a more complicated manner. Although accurate numerical treatments (Ma & Bertschinger 
1995; Dodclson et al. 1996a) exist, these complications have meant that fitting formula ( Holtzmarj 1989; 



Klypin et al. 1993; Pogosyan & Starobinsky 1995; Ma 1996) for the power spectra of such cosmologies have 



been restricted to certain regions of parameter space, e.g. fixed baryon content and critical density overall. 



In Hu & Eisenstein (1997, hereafter HE97), we showed that the late-time evolution of perturbations 



in a mixed dark matter (MDM) cosmology with CDM, baryons, and massive neutrinos could be accurately 
treated using a scale-dependent growth function. The transfer function then becomes the product of this 
growth function and a time-independent master function that represents the perturbations around recombi- 
nation. Moreover, the small-scale limit of this master function can be calculated analytically as a function 



of cosmological parameters ( HE97 ) 



In this paper, we give an accurate fitting formula for the master function. This in turn produces a 
fitting formula for the transfer functions of adiabatic cosmologies as functions of matter density, baryon and 
neutrino fractions, cosmological constant, Hubble constant, redshift, and the number of degenerate massive 
neutrino species. For the central region of parameter space, i.e. only moderate deviations from the pure- 
CDM model, the formula is accurate to better than 5% in the transfer function (10% in power). The formula 
does not attempt to fit the acoustic oscillations created by large baryon fractions but provides a good match 
to the underlying smooth function. Hence, the formula loses accuracy for baryon fractions exceeding 30%. 
Similarly, the formula has larger errors for cosmologies with massive neutrino fractions exceeding 30% or 
with matter densities outside the range 0.06 <; iloh 2 <J 0.4. In this range, however, if the baryon fraction is 
less than 10%, the accuracy improves to better than 3%. 

The outline of the paper is as follows. In §||, we review the basic results of linear perturbation theory. 
We then present the fitting function in §|| and a user's guide in As illustrative examples of the utility of 
the formula in exploring parameter space, we examine constraints on large-scale structure and high-redshift 
objects in §0. We conclude in §|[ 



2. Description of Physical Situation 

We consider linear adiabatic perturbations around a Friedmann-Robertson- Walker metric for cosmolo- 
gies with several species of particles: photons, baryons (i.e. nuclei and electrons), massive and massless 
neutrinos, and cold dark matter. The interaction between this diverse set of particles can lead to complex 
phenomenology in the growth of perturbations even in linear theory (e.g. Peebles) 1993). 



Nevertheless, the underlying physical situation remains simple. Perturbations under the so-called Jeans 
scale are not subject to gravitational instability due to pressure support or, in the case of collisionless 
particles, sufficient rms velocity. Above this scale, perturbations grow at the same rate regardless of scale. 
In general, the Jeans scale of each gravitating species is imprinted into the power spectrum. It is conventional 
to define the transfer function as the ratio of time-integrated growth on a particular scale as compared to 
that on scales far larger than the Jeans scale. 



For a relativistic species, the Jeans scale grows with the particle horizon. In a universe with cold dark 
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matter and radiation only, the Jeans scale of the total system grows to the size of the horizon at matter- 
radiation equality and then shrinks to zero as the universe becomes matter-dominated. The result is that 
the transfer function turns over at the scale of the horizon at equality. Moreover, well after equality, the 
Jeans scale has dropped sufficiently that all scales of interest are above it and hence grow at the same rate. 
The familiar result is that the spectrum of fluctuations can be written at low redshifts as a scale-independent 
growth factor times a function of scale that depends only on the size of the horizon at equality. 

With the inclusion of the baryons, another scale is imprinted in the transfer function. The baryons 
are dynamically coupled to the photons until the end of the Compton drag epoch, close to recombination 
for the standard thermal history. Prior to this time, the baryonic Jeans scale tracks the horizon (or more 
properly the sound horizon, accounting for the fact that baryons contribute mass but little pressure). After 
recombination, the Jeans scale of the baryons drops precipitously to scales smaller than those of interest for 



large-scale structure (for sub- Jeans perturbations, see Yamamoto et al. 1997). The sound horizon at the end 



of the Compton drag epoch is thereby imprinted in the transfer function in the form of Jeans or acoustic 



oscillations (c.f. EH97 ). Again, perturbations at low redshifts grow at the same rate on all relevant scales. 



The same principles may be applied to massive neutrinos. At sufficiently high temperatures, even mas- 
sive neutrinos behave as radiation; therefore their Jeans scale grows with the horizon. As their temperature 
drops with the expansion, they become non-relativistic and their Jeans scale shrinks. Physically, their mo- 
menta decay with the expansion and eventually become small enough to allow them to cluster with the 
non-relativistic matter. For this reason, the Jeans scale is often called the free-streaming scale. By the same 
arguments as before, the maximal free-streaming scale is imprinted in the transfer function. What makes 
the massive neutrino case more complicated than the cold dark matter and radiation case is that for eV 
mass neutrinos the free-streaming scale today lies in the regime of large-scale structure measurements. Thus 
the growth of fluctuations in the regime of interest is not independent of scale even at low redshifts. 

Let us examine the growth in more detail. A given scale begins outside the free-streaming length; here 
the neutrino density perturbation traces those of the other species. If the scale is below the maximal free- 
streaming scale, it will at some point cross the free-streaming scale. While in the free-streaming regime, 
perturbations in the neutrinos damp out collisionlessly while those in the cold dark matter and baryons grow 
more slowly due to the loss of a gravitating source. As the neutrinos slow down and their Jeans scale shrinks, 
the scale in question eventually crosses back out of the free-streaming regime. At this time, the neutrinos 
fall back into the potential wells of the other species and the growth rate is boosted back to its original rate. 
Even at low redshifts, some scales are still in the free-streaming regime; hence, the temporal and spatial 
dependence of the transfer function cannot be separated as before. 

If all of the massive neutrinos had the same momentum, then one could hope to describe the free- 
streaming situation more exactly, but of course the neutrinos have a thermal distribution, frozen in when the 
universe had a temperature of about 1 MeV. Hence, the transition between free-streaming and infall occurs 
smoothly and requires a Boltzmann code to follow ( Ma fc Bertschinger] 1995). HE97 showed that the result 



could be well fit by a scale-dependent growth rate; we will use this here to separate the time dependence 
from the complications of the spatial dependences. 
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3. Fitting Form 

3.1. Scales and Notation 

We begin by describing our notation. The density of cold dark matter, baryons, and massive neutrinos, 
in units of the critical density, are denoted O c , 0;,, and 0„, respectively. The total matter density is then 
flo = O c + fib + Vl v . f c , ft, and f v are the ratio of the density of these species to the total Oo. We use 
multiple subscripts to indicate summation, so that, e.g., f c b = fc + ft = (O c + Oi,)/Oo- The contribution 
of a cosmological constant A is written as Oa = A/'SHq and is not included in Orj. The Hubble constant is 
parameterized as Hq — 100/ikm s^ 1 Mpc -1 . The CMB temperature is given by Tcmb = 2.702.7K; the best 



determination to date is 2.728 ± 0.004 K (Fixsen et al. 1996; 95% confidence interval), at which it is fixed 
for most of our expressions. 

We assume that there are three species of neutrinos with a temperature equal to (4/11) 1 / 3 of the photons 
while they are relativistic. One or more of the species may be sufficiently massive to influence cosmology, 
but we only study the case where the most massive species have essentially equal masses. Then N v is the 



number of these species, and the mass of each is m u = 91 .50„/i 2 N „ 1 eV/c 2 (Kolb & Turner 1990) 



We generally work with the redshift z as our time coordinate. The redshift of matter-radiation equality 



2.50 x 1O 4 O O ^ 2 0^.7- 



(1) 



The baryons are released from the Compton drag of the photons near recombination at a redshift (sec Hu 



fe Sugiyama 1996; EH97) 



z d = 1291 



62 



(O /i 2 ) c 



1 + 0.659(O /i 2 ) 



0.82s 



[1+hiflbh 2 ) 



2^21 



(2) 



0.313(O /i 2 
0.238(O fr 2 



-0.419 



l + 0.607(O /i 



2\0.674i 



However, it is more convenient to refer this quantity to the expansion factor since matter-radiation equality, 
so we define 



Vd 



1 + Zeq 

l + z d 



(3) 



The comoving distance that a sound wave can propagate prior to z d is called the sound horizon and is (EH97 
1997) 

44.5 01(9.83/00^ 



Mpc 



(4) 



(note that the units are Mpc, not h 1 Mpc). 



As we are using linear perturbation theory, it is appropriate to work in Fourier space, where the transfer 
function depends on the comoving wavenumber k. We often parameterize k relative to the scale that crosses 
the horizon at matter-radiation equality, so as to define 



q = 



Mpc 



-oi 7 (n h z 



19.0 



(Qoffo 2 r 1/2 (l + *e q )- 1/2 . 



(5) 



"Although this is not well-defined in cases with f2„ ^ 0, we justify our choice in HE97 
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3.2. Free-Streaming and Infall 



As shown in HE97, one can decompose the transfer function into a scale-dependent growth function that 
incorporates all post-recombination effects and a time-independent master function that reflects conditions at 
the drag epoch. Hence, we write the transfer function of the density-weighted CDM and baryon perturbations 
as 

T cb {q, z) = T maBtCT (q)D cb (q, z)/Di(z) (6) 
and that of the density-weighted CDM, baryon, and neutrino perturbations as 



T c bv(q, z) = T masteI (q)D cbv (q, z)/Di{z). 



(7) 



Here, D\ is the growth factor for the universe in the absence of neutrino free-streaming (i.e. on very large 
scales), D c b and are the scale-dependent MDM growth functions, and T mastC r is the time-independent 
master function. We describe these now in turn. 



In the absence of free-streaming, the growth function takes on the usual form (Heath 1977; Peebles 



1980) 



Di{z) 
9 2 { Z ) 



5^o 



(1 + z cq )g(z) 



g(z') 3 



dz' , 



n (i + zf + (i - o - n A )(i + zf + n A . 



(8) 
(9) 



We have chosen the normalization to be D\ = (1 + 2 e q)/(l + z) at early times. For an £Iq = 1 universe, 
equation (H) yields D\ = (l + z cq )/(l + z) at all times; closed-form expressions are also available for universes 
without A (Weinberg 1972; [Edwards fc Heath 1976; Groth & Peebles 1975) and flat low-density universes 



with A (Bildhauer et al. 1992). Alternatively, one may use the fitting form (Lahav et al. 1991; Carroll et al 
1992) 



Di{z) 

n(z) 
n A (z) 



1 + z eq \ 5fl(z 



l + z 



{ti(zf/ 7 - n A (z) + [i + n{z)/2][i + r> A (z)/70]} 



(10) 



n (l + z) 3 g- 2 (z), 
^Ag- 2 (z), 



where g(z) is defined in equation (|9|). 

The presence of neutrinos suppresses the growth of fluctuations on scales sufficiently small that the 
neutrinos' velocity allows them to escape the perturbation. This alters the logarithmic growth rate (Bond 

(11) 



et al. 1980) according to the factor (i = cb,c) 



l r 



5 - VT+ 24/i 



> 0. 



Then the growth rates in the presence of free-streaming are ( HE97 ) 



and 



D cb (z,q) 



D chv {z,q) 



Di(z) 



Jcb 



y{ S {q;fu) 

Di(z) 



0.7 



Pcb/0.7 



D x {z 



,1-Pc 



0.7 



Pcb/0.7 



Di{z) 



l-Pc 



(12) 



(13) 
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for the CDM+baryon and CDM+baryon+neutrino cases, respectively. In both cases, the free-streaming 
epoch as a function of scale is 

y is (q) = 17.2/„(1 + 0.488/- 7 / 6 )(iW/„) 2 . (14) 

Note that increasing N v at fixed £l u prolongs free-streaming by making the neutrinos less massive and hence 
faster moving. 

The functions D c b and D c b u contain all the dependence of the transfer functions on time, curvature, 
and cosmological constant, and moreover relate the two transfer functions to a single master function. Of 
course, in a cosmology with no massive neutrinos, D c b = D c b u = D\ and the master function is simply the 
usual post-recombination transfer function. 



3.3. The Master Function 

The master function reflects the spectrum of perturbations at the drag epoch. As such, it can only 
depend on VLqU 2 : fb, fv, and N v . In the case without massive neutrinos, this reduces to the CDM+baryon 



transfer function. EH97 described the phenomenology of this function. In particular, the presence of baryons 



suppresses power on scales smaller than the sound horizon at the drag epoch and introduces a series of 
oscillations in the transfer function that damp as one moves to smaller scales. For moderate baryon fractions 
(/b <; flgh 2 + 0.2), the oscillations are fairly small but the suppression can be important (roughly 5/& in 
power). Because of the increased complexity of adding massive neutrinos to the form, we opt not to fit 
the oscillations and instead tailor a formula that runs through the center of the oscillations and properly 
incorporates the small-scale suppression. Of course, this means that the formula will not be appropriate for 
cases where the oscillations are large, roughly f2;,/(f2 c + Of,) ;> floh 2 + 0.2. 

On scales much larger than the sound horizon at the drag epoch and the horizon at the time when the 
neutrinos finally become non-relativistic, the effects of pressure fluctuations and collisionless damping are 
not important. Therefore, the transfer function will match that of a pure-CDM model. On small scales, 
we have solved the evolution equations analytically and therefore can calculate the amount of small-scale 



suppression due to the baryons and neutrinos (HE97). We use these two limits to anchor our fitting form. 



The suppression of power in the master function on small-scales is primarily due to baryons, although 
neutrinos do contribute a residual coefficient not included in the free-streaming growth function 3 . The 



amount of small-scale suppression is given as (HE97) 



<xv(fu,fb,y<i) = 



f c 5 - 2(p c +p cb ) 



1 - 0.553/„ fc + 0.126/, 



f, 



cb 



1 + 



4p c b 

- Pcb 



1 



1 - 0.l93y%N; + O.I69/^o- 2 



(15) 



2 V' ' (3-4p c )(7-4p cb ), 
We choose to include this suppression by a scale-dependent rescaling of the zero-baryon shape parameter T 



(c.f. EH97). The suppression occurs rapidly near the sound horizon (defined in eq. H): 



- off 



= fink 2 



1 - 



1 + (0.43fcs) 4 



(16) 



s The total suppression of T c j on small-scales is a v D 1 



Then we use this effective wavcnumber in the zero-baryon form, 



W*) = j -Jl^, (18) 

L = ln(e+ IMpcy/c^qes), (19) 
325 

° = 144 + 1 + 60.5^ ' (20) 

& = (1 - 0.949/U)- 1 , (21) 

to produce a form that breaks from the large-scale, pure-CDM formula to the small-scale solution. 

We find, however, that this formula is inaccurate around the scale of the horizon at the epoch when 
the neutrinos slowed to non-relativistic speeds, the so-called maximal free streaming scale (cf. This is 
because the form of yt s (eq. Q) assumes that the neutrino velocity scales simply as v cx (1 + z). In fact, 
it is the momentum that carries this scaling, while the velocity cannot exceed c. This error causes us to 
overestimate the maximal free-streaming scale (note that in eq. juj the running of this scale with redshift is 
cutoff at the equality epoch). In turn, the growth functions D c b and D c \, v provide too much free-streaming 
suppression on these scales, although since the scales are well above the free-streaming scale for z <; 30, no 
spurious time-dependence is introduced at late times. 

For /„ < 0.3, this error can be fixed by the following multiplicative correction: 

B(k) = 1+ - ^.e " , (22) 



ft 




k 

q v = == — = 3.92?W^. (23) 



The master function is then 

T mastcr (fc) = T sup {k)B(k). (24) 



3.4. Performance 

In Figures 1 and 2, we compare the fitting formula to the numerical evaluation (using the CMBfast 



code of Seljak & Zaldarriaga 1996 v. 2.3) of the transfer function. Figure 1 shows the most common MDM 
model — f2o = 1, h = 0.5, £!& = 0.05, Vl u = 0.2, and N v = 1 — at two different redshifts. Figure 2 shows other 
cases — high baryon fraction, low Slo, N v — 2, and high iloh 2 — at redshift zero. 

For 0.06 £ VL Q h 2 £ 0.40, n b /Q, < 0.3, fV^o < 0.3, z = 0, and N v = 1, the accuracy of the fitting- 
formula is quite high. For baryon fractions of 5%, the acoustic oscillations are small and the fit is better than 
2% on all scales. At higher baryon fractions, the oscillations become more prominent, and so the maximum 
level of the residuals grows, although the residuals would at least partially cancel for many applications. 
Performance for /t, < 0.3 is nearly always better than 5% (and often <J 3%) when comparing to the non- 
oscillatory portion of the transfer function. The small-scale fit for q > 0.25 (k > 1 Mpc -1 ) is always better 
than 2% accurate. Behavior for N v — 2 is similar, although different numerical codes seem to be inconsistent 
for q 3> 1 at the several percent level. Performance at z = 9 is at most 1% worse than that at z = 0; at 



- 8 - 




Fig. 1. — Comparison of the fitting formula to the numerical results of CMBfast ( Seljak fc Zaldarriaga 



1996). (a) Results at z = 0. (b) z — 9. Upper panels: Transfer functions divided by a fiducial purc- 



CDM transfer function, formed by using equation (18) with a„ = (5 C = 1 and g c ff = q. Density-weighted 



CDM+baryon (short- dashed) and CDM+baryon+ncutrino (long-dashed) transfer functions arc shown. Lower 
panels: Fractional residuals. The cosmology is flo = 1, fi„ = 0.2, Qb = 0.05, and h — 0.5. 

z = 29, performance can degrade by 4% at the lowest values of iloh 2 (where z/z cq is its largest). Hence, 5% 
accuracy is achieved only for z < 30. 

For fl^h 2 ;> 0.4, the fitting formula tends to overestimate the transfer function on scales just below 
the sound horizon (k rss 0.3 Mpc -1 ) and underestimate it just above the sound horizon (k 0.1 Mpc -1 ); 
we display this problem in Figure ||<i These errors are only a few percent for /& < 0.1, but grow to 10% 
by fb ~ 0.3. The culprit is our reliance on using an effective T within a pure-CDM transfer function; for 
high floh 2 the sound horizon corresponds to q <C 1, thereby altering the portion of the curve we use for our 
transition (eq. |]l6| ). For flgh 2 <J 0.06, the opposite situation occurs; moreover, the power series in 1 + yd in 
equation ( |l5| ) becomes less accurate. For fl h 2 = 0.025, the errors are 5% for both low and moderate baryon 
fractions. 

For models with fli,/(flb + f2 c ) ^ ^o^ 2 + 0.2, the baryon oscillations exceed 10% in amplitude, which 
may be a problem for some applications. By floh 2 + 0.4, the oscillations are of order 40% (c.f. Fig. 5 of 



EH97). We note that the location of the peaks in wavenumber seems essentially constant with varying i\ 



thus the formulae in EH97 could be used to give the location but not the amplitude. 



For neutrino fractions exceeding 30%, our correction for the behavior near the maximal free-streaming 
scale (eq. p2]| ) is too small, leading to significant errors (8% for f v sa 0.5, increasing thereafter). 

For ilo < 1) we have only tested the fit on intermediate scales for flat universes (i.e. £7a = 1 — Slo)- We 
have tested the small-scale limit in both open and flat cases and found excellent accuracy. 

Note that our formula works equally well for cases with f2„ = 0. Of course, in this case, should one wish 



to fit the baryon oscillations, one should use the fitting formula in EH97 
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Fig. 2. — Same as Figure |l], but for different choices of cosmology, (a) High-baryon model, with SI5 = 0.3, 
fi„ = 0.2, fi = 1, and h = 0.5. (b) Low-density, flat model, with Q = 0.3, Q A = 0.7, fi„ = 0.06, Q b = 0.03, 
and h = 0.5. (c) Model with two degenerate neutrino species (N v = 2). f2& = 0.05, il„ = 0.2, ilo — 1, and 
h = 0.5. (d) Model with high n h 2 , beyond our range of [0.06,0.40]. n b = 0.2, VL V = 0.2, n Q = 1, and 
h = 0.8. 



4. User's Guide 



We present here a user's guide to the fitting formulae of the previous section. The fitting formula for 
the density- weighted matter transfer function, with and without neutrinos, is given by equations (|l]) ( |24|) . 
For cases with Vl v 7^ 0, these functions are time-dependent and involve the growth factors in equations (|l2j) 
and (|T5|). We now detail how to use the transfer function to construct power spectra and measures of mass 
fluctuations. 
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4.1. Power Spectra 

The power spectrum is constructed from the transfer functions in the usual way: 

k 3 / rk\ 3+ " 

—sP(k, z) = 5 H (—) T 2 (fc, z)Dl(z)/Dl{Q). (25) 



2tt 2 v,j u \H 0/ 

where 5h is the amplitude of perturbations on the horizon scale today, and n is the initial power spectrum 
index, equal to 1 for a scale-invariant spectrum. Note that the usual growth function D\ from equation (||) 
is used, not D cb or D cbv . 

In cases with Q u ^ 0, there are three transfer functions and hence three power spectra that may be 
constructed. Using T cb from equation (^) in equation ( p5|) yields P cb , the power spectra for the CDM and 
baryons. Likewise using T cbv from equation (Q) yields P c bv, the density- weighted power spectrum of the 
CDM, baryons, and massive neutrinos. The power spectrum of the massive neutrinos themselves can be 
constructed from the functions above. One can subtract the transfer functions given above to find the 
transfer function for the neutrinos alone: 

T v = f-\T cbv - f cb T cb ). (26) 

On small scales, this function goes to zero, but we have not carefully modeled this. Thus, T v has density- 
weighted errors similar to those in T cb ; that is, the residuals in the fit for f v T v will be similar to those of T cb . 
This transfer function may be employed in equation ( |25| ) to obtain P V) the power spectrum of the massive 
neutrinos. 

Power spectra for velocity fields can be similarly obtained by considering the continuity equation, which 
relates them to time-derivatives of the density fluctuations. It is standard to express this in terms of the 
quantity / = —d\og(D)/d\og(l + z), where D is the growth function. In a model with massive neutrinos, 
this becomes a scale-dependent quantity. We can differentiate equation ([l2]) directly to find 

fcbik, z) = f (z) J 1 ^ I (27) 

with fo(z) as the value of / in the absence of free-streaming: 



s -SEgTTT) ^ n < 2) "'' + m Q -^ 1 1 + ^ I ■ < 28 » 



where the approximation (Lahav et al. 1991) uses Q(z) and Q\(z) from equation ([H 



The power spectrum for the velocity field for the CDM and baryons is then 

2 



^M)=f /c t' z) y r i p*v<,*), (29) 



(1 + z)k 

where g(z) was defined in equation (O). A similar relation follows for the velocity field of the density- weighted 



matter. 



4.2. COBE Normalization 

To normalize the power spectrum to the COBE-DMR measurement, one may use the fitting formulae 



of Bunn fc White (1997) to fix 5h- For cases with no CMB anisotropics from gravitational waves, one has 

S H = 1.94 X l -5 r , o 0.785-0.051„flo e -0.95n-0.169n^ A = 1 - Q , (30) 
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6 



H 



1.95 x 



5o-0-35-0.191nfio-0.17S -fi-0.14n 2 



(I 



A = 0, 



(31) 



valid for 0.7 < n < 1.2. For flat cosmologies with the gravitational wave contributions of power-law inflation 
(which requires n < 1), 5h is 

S H = 1.94 X 10-5ft-0.785-0.051nfio e n+1.97n 2 i A = 1 - O . (32) 

For open cosmologies with the minimal gravitational wave contribution from power-law inflation (again, 
n < 1), one has flHu fc White] 1997) 



5 H = 1.95 X 10-8 n -0.86-0.191nno-0.16ft e -1.02ft-1.70ft ^ A = Q (33) 

In all cases, n = n — 1 and the fits extend from 0.2 < fio < 1. The la statistical uncertainty in the 
COBE-normalization is 7%, primarily due to cosmic variance. 



4.3. Mass Fluctuation Measures 

The rms amplitude of mass fluctuations inside a particular spherically-symmetric window is 

dk k 3 



k 2n^ k) 



1/2 



(34) 



where P(k) is the power spectrum and Wn{k) is the Fourier transform of the real-space window function 
Wfi(r). Either P c i, or may be used depending on the application. The two most popular choices for 
window functions are the real-space spherical tophat of radius R: 



W R (r) oc 
W R (k) = 

Mr = 



1 if r < R, 
otherwise, 
3 



(fci?) 3 
^ Pc n R\ 



(sin kR — kR cos kR) , 



and the Gaussian window of scale length R: 

Wfl(r) oc exp 



W R (k) 



exp 



r 2 \ 
2R?) ' 
(kR) 2 



M R = (2ir)V 2 p c n R 3 . 
Here, M R is the mass included in the window. 



(35) 
(36) 
(37) 

(38) 

(39) 
(40) 



5. Observational Constraints 



The fitting formula presented in §|3| allows one to manipulate statistics of the power spectrum as functions 
of cosmological parameters much more easily than a suite of Boltzmann integrations would allow. As 
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examples of its use, we consider predictions for the power spectrum of large-scale structure, the abundance 
of clusters of galaxies, damped Lya systems, and the Lya forest. The theoretical power spectrum is related 



to these observations via the rms amplitude of mass fluctuations (§4.3) on various scales. 



We consider 2-dimensional cross-sections in parameter space by varying the baryon and massive neutrino 
fractions in four fiducial models: the standard case of Slo = 1, h = 0.5, n = 1, and N v = 1; a tilted variant 
with n = 0.95 and a tensor contribution to the CMB anisotropy; a variant with a second neutrino species 
(N v = 2); and a low-density flat universe with fig = 0.35, h = 0.7, n — 1, and N v — 1. We choose these 
models in order to explore the various ways of addressing what has been identified as the key problem of 
standard CDM (i.e. Qo = 1, h = 0.5, n = 1, and trace or zero baryon and neutrino content), namely the 
overproduction of power on galaxy and cluster scales relative to larger scales (e.g. |Efstathiou ct al. 1992; 
Ostriker 1993; Dodclson ct al 1996b). As explained i n §[| adding massive n eutrinos (Bchaefer & Shafi 1992; 
Davis et al, 1992; Taylor & Rowan-Robinson 1992; Holtzman fc Primack] 1992) or baryons (White et al 



1996) reduces small-scale power. This alone may be sufficient to satisfy constraints. However, other simple 
extensions act to suppress power and may produce a better fit to the data. Adding a red tilt (n < 1) to the 
initial powe r spectrum (see e.g. | Ccn ct al. 1992) or lowering the density parameter Qq (see e.g. jEfstathiou 
at al. 1992; Ostriker & Steinhardt 1995) are common approaches. Here, we also consider the addition of a 



second species of massive neutrinos (Primack ct al. 1995); this helps because it further reduces power on 
cluster scales while leaving the small-scale power essentially unchanged. 



5.1. Power Spectrum Shape 



We begin at large scales and consider the shape of the linear power spectrum as reconstructed from 



galaxy surveys. Peacock fc Dodds (1994) considered a collection of data sets vis-a-vis zero-baryon power 
spectrum models; they found that scale-invariant models with T = ttoh = 0.255 ± 0.017 provided the best 
fit. However, adding baryons and/or neutrinos alters the shape and hence the best fit. We do not perform 
this detailed reanalysis here; rather we use the ratio of large-scale to small-scale power as a proxy for the 



shape (see e.g. White et al. 1996). In particular, we construct the ratio of the amplitude of density- weighted 
fluctuations (P c bu) within a 50/i _1 Mpc tophat to those within a 8/i _1 Mpc tophat, i.e. crso/cs- The range 
r = 0.25 ± 0.05, which we conservatively adopt, converts to crso/crs — 0.151 ± 0.016. Note that high values 
of r produce lower values of a^o/vs- 

We display this ratio in Figure [|a-d (left panel) as a function of baryon fraction and neutrino fraction. It 
is important to note that baryons play as important a role as neutrinos in suppressing power on cluster scales 



relative to larger scales ( White et al. 1996). For example, even the cosmic concordance model ( Ostriker fe 
Steinhardt 1995) of Qq = 0.35 with h = 0.7, which seems to have an appropriate T, stretches the constraint 



when pushed to the 10-15% baryon fraction suggested by cluster mass determinations (e.g. |White ct al 
1993b; |David et al.| 1995; |White fe Fabiarj 1995; |Evrard| 1997). 



5.2. Cluster Abundance 



The present-day abundance of rich clusters of galaxies is a sensitive probe of mass fluctuations on the 
8/1- 1 Mpc scale (|Evrard| 1989; |White et al.| 1993a; |Ekc ct al.| 1996; |Viana fc Lid"dTej 1996; |Bond fc Myers) 1996; 

(1997) 

0.65 



Pen 1997). We adopt the determination of 



CTg PS 0.5iln 



(41) 
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Fig. 3. — (a-d) Ratio of fluctuations within 50h 1 Mpc spheres to that within 8h 1 Mpc spheres as a function 
of neutrino and baryon fractions for various cosmologies. The shaded region shows the preferred range of 
<Wcr 8 = 0.151 ± 0.016. (a) Q = 1, h = 0.5, n = 1, N v = 1. (6) As (a), but with n = 0.95 and tensors 
as per power-law inflation, (c) As (a), but with two degenerate neutrino species (N u = 2). (d) f2o = 0.35, 
57a = 0.65, h = 0.7, n = 1, N v = 1. All are COBE-normalized. (e-h) Amplitude of fluctuations within 
8/i _1 Mpc spheres for the cosmologies given in (a-d), respectively. The shaded region is the preferred range 
of a s 



(0.5 ± 0.15)1V - 65 ( Pen 1997) and a$ > 0.59 (Fan et al. 1997). 



and take a conservative range of 30% errors (i.e. 0.5 ± 0.15). 

The time evolution of the abundance of clusters provides a way to isolate as from its f^o dependence. By 



comparing the abundance of high-redshift clusters relative to the present-day abundance, Fan et al. (1997) 
found as = 0.83 ± 0.15(ler). To be conservative, we employ a 2a lower limit (from the relevant quantity, 
erg" 2 ) of a s > 0.59. 

Fig. ||e-/i (right panel) shows the cluster abundance constraints for the same 4 models as Fig. ||a-d. 
As is well-known, high-f^o cosmologies with small tilt and trace baryon and neutrino content overproduce 
present-day clusters. Adding a substantial fraction of baryons or neutrinos makes the models marginally 
consistent with both the present-day and high-redshift cluster abundances. 



5.3. Damped Lya Systems 



Cosmologies with moderate neutrino fractions have a strong suppression of power on small scales. This 
implies that they form proto-galactic systems later than pure-CDM models. Indeed, these models may have 
trouble forming high-redshift objects such as quasars (e.g. Ma & Bertschingei 1994; Liddle et al. 1996), 
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Fig. 4. — (a-d) Amplitude of fluctuations at z = 4 within a Gaussian window of mass corresponding to 
halos of 50 km s _1 Mpc -1 circular velocity. Cosmologies are as per Figure |[ The shaded region indicates 
cosmologies where the neutral gas in halos of v c > 50 km s _1 Mpc~ x [using the prescription of Klypin ct al. 
(1995) and /m = 1] exceeds that observed in damped Lya systems. The region to the left of the dashed 
line is the allowed region for /hi = 0.1. (e-h) Amplitude of fluctuations at z — 3 for a Gaussian window of 
radius 0.0416(ft ^ 2 )" 1/2 Mpc, suggested by pnedinj (1997) as a indicator of the slope of the column-density 
distribution of the Lya forest. The region OLya > 2 is shaded. Cosmologies are as in Figure |[ 



UV-dropout galaxies ( Mo fc Fukugitaj 1996), and damped Lya systems ( Mo fc Miralda-Escude 1994; Ma fc 



Bcrtschinger 1994; Kauffmann & Chariot 1994; Klypin et al. 1995). We focus on the last of these. 



Observations of damped Lya absorption systems in QSO spectra may be interpreted as a measurement 
of the mean density of neutral hydrogen f2 gas , in units of the critical density. Recent measurements at z = 4 
( [Btorric-Lombardi et al. 1996) find this to be 



f\as(2 ~ 4 ) = ( 9 - 3 ± 3 - 8 ) x KT 4 ^ 1 (1 + z) 3/2 g(z) 



Assuming Poisson statistics, we adopt a 95% lower limit of 43% of the central value. 



(42) 



One can estimate an upper limit to the value of f2 gas in a particular cosmology by assuming that all gas in 



proto-galactic halos is neutral and using the Press-Schechter formalism ( Press fc Schechtei 1974) to estimate 



the number of such halos (Mo & Miralda-Escude 1994; Kauffmann & Chariot 1994; Ma & Bcrtschinger 



1994; Klypin et al. 1995; Liddle et al. 1996). These works differ in their Press-Schechter implementation; 
here we adopt the conservative assumptions of Klypin et al. (1995). We define ctdla to be the amplitude 



of fluctuations (using P c bv) inside a Gaussian window of a scale corresponding to a circular velocity v c of 
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50 km s 1 Mpc 1 . The relation between mass and velocity is ( Narayan fc White] 1987) 



M = 



89GH g(z) ' 



(43) 



where g(z) is defined in equation (^|). Then the density of neutral gas arising from all halos with velocities 
greater than 50 km s _1 Mpc -1 is 



a, 



/mSlfcerfe 



V2a_ 



DLA 



/Hi < 1 



(44) 



where /m is the fraction of neutral gas, erfc(cc) is the complimentary error function, and we take a density 
threshold of S c = 1.33. 

In Figure Qa-d (left panel), we plot odla as a function of cosmological parameters. We superpose 
the constraint implied by comparing equation ( fi"2|) to equation ( f44|) with /m = 1. As found by previous 
studies, MDM models with /„ ;> 0.3 underproduce high-redshift halos; the constraints are tighter for higher 
fb, red-tilted, and degenerate-neutrino models. The limits in Figure^ are actually extremely conservative; 
hydrodynamical studies (Ma et al. 1997; Gardner et al. 1997a, b) infer /m ^ 0.1 in tested cases. Correspond- 
ingly, we plot the limits if /hi = 0.1 in dashed lines to show how this uncertainty affects the cosmological 
constraints. Since it is difficult to scale these numerical corrections as functions of cosmological parameters, 
especially if varying Oj, (Gardner et al. 1997b), the /hi = 0.1 line should not be taken as a firm constraint. 



5.4. Lya Forest 



If the low column density absorption features in QSO spectra arise from mild density and velocity 



perturbatio ns in the IGM flCen et ajj 1994; [Petijean et al| 1995; [Zhang et al] 1995, 1997; [Miralda-Escudd 
ct al. 1996; Hernquist et al. 1996), then the correlations and column-density distribution of the lines may 



yield robust information about the power spectrum on sub-Mpc scales. Croft et al. (1997) demonstrated 
that the power spectrum of simulations could be reconstructed from absorption spectra drawn from them. 



Gnedin (1997) showed that the power-law exponent of the column density distribution in various cosmological 
simulations is strongly correlated with the amplitude of linear fluctuations on the smallest collapsing scales. 
Comparing to the observed distribution suggests a lower bound on the amplitude of fluctuations on mass 
scales near 10 9 Mq at z = 3. In particular, the quantity CLya, defined as the fluctuations inside a Gaussian 
window of radius R = 0.0416(rio^ 2 ) -1 ^ 2 Mpc using P c b, is constrained to be greater than 2.0 at z = 3. The 
scale is chosen to approximate the Jeans length at z = 3 for common thermal histories ( pnedin| 1997). This 
constraint is plotted in Figure |l|e-/i (right panel) . 



5.5. Summary 

Even in the standard COBE-normalizcd, n = 1, Qq = 1 model, the inclusion of a moderate fraction of 
baryons or neutrinos can decrease as to an appropriate level. Models that accomplish this by the neutrino 
fraction alone produce insufficient power to explain damped Lya absorption systems. Models that accomplish 
this by the baryon fraction alone require baryon densities far in excess of big bang nucleosynthesis predictions. 
Although a compromise of £1/, = 0.15 and f2„ = 0.25 would work, Figure|| shows that no model in this scenario 
fits the <75o /<J8 an d <y& constraints. 
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However, the modest change of altering the tilt to 0.95 (with tensors) or adding a second degenerate 
neutrino species opens regions of parameter space that match both large-scale structure and constraints from 
damped Lya systems. For example, models with fij = 0.1, h = 0.5, and either f2„ = 0.15 with n — 0.95 or 
il v = 0.2 with N v = 2 produce <7 8 « 0.64 and match the high-redshift constraints with /hi ~ 0.2. /hi ~ 0.1 
could be achieved by reducing the COBE normalization by 7% (a 1 — a shift) and decreasing f2„ so as to keep 
<j% constant. The higher value of — as compared to the canonical value of 0.05 from Walker et al. (1991) 
but in agreement with Tytler et al. (1996) — is doubly useful in meeting the requirements: the suppression 



due to baryons occurs at larger scales than that from neutrinos and therefore alters a% more effectively, while 
the additional baryons are available to produce high-redshift absorption. Moreover, this value of fib better 



agrees with the baryon fraction in clusters (e.g. White et al. 1993b; David et al 



Evrard| 1997) and that inferred from the Lya forest (e.g. [Miralda-Escude et al 



1995; |White fc Fabian 1995; 
1996; IWeinberg et al. 1997: 



Zhang et al. 1997) 



The more common solution to the problems of f^o = 1 CDM is to reduce the value of £7 to around 0.3. 
As shown in Figures [3] and |4| this satisfies the quoted constraints when used with small baryon and neutrino 
fractions. An important lesson of the figures, however, is that small admixtures of baryons or neutrinos 
can make significant changes. For example, the Qq = 0.35 flat model presented here tends to underproduce 
power even with only a 10% baryon fraction (n^h 2 = 0.017). The lack of small-scale power in such models 
places more stringent limits on f2„/f2o than in high-density cosmologies; this implies a stronger limit for the 
neutrino mass ra w oc Q„h 2 . Of course, a small blue tilt would help the situation. 



It is very interesting to note that early results from modeling the Lya forest (§5.4, Fig. ^e-/i) ar e more 
effective at excluding models than constraints from damped systems. The limits suggested by Gncdin (1997) 
would eliminate all of the Slo = 1 models studied in Figures || and|i[ Perhaps models with blue tilts (n > 1) 
would succeed in producing sufficient amounts of small-scale power, although of course they would require 
more suppression of power at cluster scales relative to COBE. Constraints from the Lya forest are still only 
preliminary, but they appear to be quite promising. 



6. Conclusion 

In this paper, we have considered adiabatic models composed of baryons, cold dark matter, and massive 
neutrinos. We have presented a fitting formula for the linear transfer function of such models, including the 
possibility of fio 7^ 1 and multiple degenerate neutrino models. The parameter space covered by the formula 
is much larger than that previously available; we provide functions of space, time, and six cosmological 
parameters. The accuracy is <J 5% in the central range of 0.06 <J £loh 2 <; 0.40, Slfc/f^o < 0.3, VL v /Q,q < 0.3, 
and z < 30 and improves to < 3% for models with baryon fractions below 10%. 

An accurate, general fitting formula allows one to calculate statistics of the power spectrum as func- 
tions of cosmological parameters quite efficiently. As an example of this, we presented several different 
observational tests and displayed the constraints as functions of baryon fraction and neutrino fraction for 
various choices of the other cosmological parameters. Baryons and neutrinos are both effective at suppressing 
small-scale power relative to that on larger scales. We find that models with large baryon fractions are less 
"observationally-challenged" , in that a given reduction on cluster scales (i.e. erg) imposes less suppression on 
very small scales, where power is needed to produce damped Lya systems and other high-redshift objects. 



Hence, if Qj, is as large as 0.1, as suggested by Tytler et al. (1996) for h = 0.5, then constraints on mixed 



dark matter models are weakened. For example, with a tilt of n = 0.95, an MDM model with Sl„ = 0.15 
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and fib — 0.1 fares significantly better than one with fl v = 0.2 and fib — 0.05. Finally, we find that the 
constraint on the small-scale power as derived from the slope of the column-density distribution of the Lya 
forest flGnedinj 1997) is an extremely powerful limit on MDM models. Further work is needed to test the 
robustness of this inference. 

The observations discussed above place constraints upon the neutrino mass, although these limits vary 
with other presently unknown parameters, e.g. flo, h, fib, and n. Future CMB observations should precisely 



determine these quantities ( Jungman et alj 1996; Bond ct al. 1997; |Zaldarriaga et al] 1997) but will have 



little leverage on f2„ ( Ma fc Bertschinger 1995; Dodelson et al. 1996a). However, the combination of CMB 
data with large-scale structure observations will allow a robust determination of fl v . Further observations 
and modeling of damped Lya systems and the Lya forest will corroborate this but may not be clean enough 
to yield a precise measurement. If flo is found to be low, our sensitivity to the neutrino mass will be stronger 
because the suppression of small-scale power depends on fl^/fl ; this differs from the trend in the CMB, 
where lowering fl shifts the effects of neutrinos to smaller angular scales. This illustrates the power of 
combining cosmological data sets with regard to determining the properties of the dark matter. 

The formulae in §0 of this paper are available in electronic form at 



http: / / www.sns.ias.edu / ^whu/transfcr / transfer, htm] . 



We also include a driver that calculates COBE-normalized erg and other constraints from 
cosmological input parameters. 



as a function of 
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